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^ I In this paper we estimate the influence of the shadowing effect and initial state parton 

Tjj- ' energy loss in the quarkonium production at HERA-B. We analyze the xp behavior of the 

^ ■ effective exponent a^xp) and present a comparison with the preliminary HERA-B data for 

^ . J /"if production. Moreover, we estimate the magnitude of these effects in the J/^ production 

> : at RHIC. 

OO , 

in ■ 

if^ I The relativistic collider facilities RHIC and LHC provide the opportunity to systematically 

^ ■ study the physics of hot and ultradense matter in hadron-nucleus (pA) and nucleus-nucleus 
Q ■ [AB) collisions at high energies (For a review see, e.g. [1, 2]). The systematic study of pA 
"^H ! collisions at the same energies is essential to gain insight into the structure of the dense medium 
P-i! effects. Such effects, as the energy loss and shadowing, are absent or small in pp collisions, but 
become increasingly prominent in pA collisions, and are of major importance in AA reactions. 
By comparing pA and AA reactions involving very heavy nuclei, one may be able to distin- 
guish basic hadronic effects that dominate the dynamics in pA collisions, from a quark-gluon 
• 1— I , formation predicted to occur in heavy ion AA collisions. To gain insight into the underlying 
^ i hadronic processes, one has to study collisions that are expected to not lead to a QGP forma- 
d I tion. Once the physics of "QCD at high densities" is better understood, the mechanisms of 
quark-gluon plasma formation and related collective phenomena in heavy ion collisions could 
be disentangled from the basic hadronic effects. 

In this paper we study the influence of the nuclear medium effects in the quarkonium 
production, particularly of the J/ip, which is one of the proposed signatures of the QCD phase 
transition [3]. In particular, we will consider the shadowing effects in the parton distributions 
and the initial state parton energy loss, which have the strongest influence on the xp {= X1—X2) 
behavior of the cross section (For a similar analyzes see, e.g. [4]). Currently, the A dependence 
of J/\E' production at > is known to rather high precision at several different energies 
(See e.g. [5]). On the other hand, the behavior of the cross section and the magnitude of 
the nuclear medium effects at negative xp region is still an open question. This situation 
should be improved by the experimental analyzes of quarkonium production in the flxed-target 
pA experiment HERA-B at DESY, which measures the quarkonium A dependence over -0.5 
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< a;^ < 0.3. First results for the J/^ and T total cross sections have been recently published 
[6] . It motivates a detailed study of the A dependence for quarkonium production considering 
the current models for the nuclear medium effects. Here we focus our analyzes in the HERA- 
B kinematical range considering two estimates for the magnitude of the initial state parton 
energy loss and two distinct parameterization for the shadowing effects. Moreover, we compare 
these predictions with the preliminary data recently obtained for the quarkonium production 
at negative values of xp [7]. As we will show, these data could be used to discriminate between 
the different models for the nuclear medium effects. As an extra possibihty of discriminate the 
different effects, we also present the corresponding predictions for RHIC energies. 

Lets start presenting a brief review about the nuclear medium effects. One of the nuclear 
medium effects is the nuclear shadowing, which is the modification of the target parton distri- 
butions so that xq'^{x, Q^) < Axq'^{x, Q^), as expected from a superposition of pp interactions 
(For a review see, e.g. [8]). In the last years several experiments have been dedicated to 
high precision measurements of deep inelastic lepton scattering (DIS) off nuclei. Experiments 
at CERN and Fermilab focus especially on the region of small values of the Bjorken variable 
X — Q'^/2Mu, where = —q'^ is the squared four-momentum transfer, u the energy transfer 
and M the nucleon mass. The data [9, 10], taken over a wide kinematic range, have shown 
that the proton and neutron structure functions are modified by a nuclear environment. The 
modifications depend on the parton momentum fraction: for momentum fractions x < 0.1 
(shadowing region) and 0.3 < a; < 0.7 (EMC region), a depletion is observed in the nuclear 
structure functions. These two regions are bridged by an enhancement known as antishadowing 
for OA < X < 0.3. We refer to the entire phenomena as the nuclear shadowing effect. 

The theoretical understanding of in the full kinematic region has progressed in recent 
years, with several models which describe the experimental data with quite success [8]. Here 
we will restrict ourselves to the descriptions which use the DGLAP evolution equations [11] to 
describe the behavior of the nuclear parton distributions. In particular, Eskola, Kolhinen and 
Salgado (EKS)[12] have shown that the experimental results [9] presenting nuclear shadowing 
effects can be described using the DGLAP evolution equations with adjusted initial parton 
distributions. The basic idea of this framework is the same as in the global analyzes of parton 
distributions in the free proton: they determine the nuclear parton densities at a wide range 
of X and > Ql = 2.25 GeV^ through their perturbative DGLAP evolution by using the 
available experimental data from lA DIS and Drell-Yan (DY) measurements in pA collisions 
as constraint. In this approach, the nuclear effects are taken into account in the initial parton 
distribution x/^(x, Ql) of the DGLAP evolution. EKS have expressed the results in terms of 
the nuclear ratios Rf{x,Q^) for each parton flavor / in a nucleus with A nucleons {A > 2), 
at 10^^ < X < 1 and 2.25 GeV"^ < < lO^GeV"^. Other groups have considered different 
set of data [13] and/or next-to-leading order corrections to the DGLAP equation and proposed 
a distinct approach [14] to describe the nuclear effects. In particular, De Florian and Sassot 
(DS) [14] proposed a framework where each nuclear parton distribution is described by a con- 
volution of the corresponding free nucleon parton distribution with a simple flavor dependent 
weight function that takes into account the nuclear effects. As pointed out by the authors, 
the Mellin transform techniques allow a straightforward parameterization of the nuclear parton 
dynamical evolution with a few parameters and an interpretation of the nuclear effects as a 



2 




Figure 1: Comparison between the EKS and DS parameterizations for the nuclear ratios 
Rll\x,Q'), <^(x,Q2) and Rl'\x,Q'). 

redistribution of longitudinal momentum among the partons in the nucleus. 

In Fig. 1 we present a comparison between the different parameterizations for the nuclear 
ratios R^^{x,Q^), R^^{x,Q^) and R^{x,Q'^) at A = 184. We can see that these parameteri- 
zations predict very distinct behavior for the nuclear parton distributions. In particular, the 
magnitude of the antishadowing in the nuclear gluon distribution is still an open question. 
This scenario should be improved by the experimental analyzes of the quarkonium production 
at HERA-B, since it probes the parton distributions in this x range. 

Another important effect in nuclear collisions is the initial state energy loss. In the last years 
the understanding of partonic energy loss has been extensively developed (For a review see e.g. 
[1, 2]), because of the expectation that the order of magnitude of the effect in hot matter is 
much larger than in cold nuclear matter, which implies that the resulting jet quenching can be 
considered a probe of the QGP formation. Our analyzes are focused in the parton energy loss in 
cold nuclear matter. It has been studied by Gavin and Milana [15] and subsequently developed 
by Brodsky and Hoyer [16] and Baier et. al. [17], considering a multiple scattering approach that 
essentially depletes the projectile parton momentum fraction, xi, as the parton moves through 
the nucleus. The basic idea is that both the quarks and gluons can scatter elastically and loose 
energy before the hard scattering. Consequently, the original projectile parton momentum 
fraction Xi when the parton first entered the target is modified to x'l = Xi — Axi, where Axi 
represents the loss in xi due to multiple scatterings, being x[ the projectile parton momentum 
fraction involved in the hard scattering. One has that the shifted value x[ enters the partonic 
cross sections but the parton distributions must be evaluated at the initial xi. Considering 
the relation between the averaged radiative energy loss —dE/dz and the characteristic squared 
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transverse momentum of the parton (p^p^), derived in Ref. [17] and given by 

dE Sag 2 



(1) 



one obtains that Axi is then 




(2) 



where La is the nuclear medium length. As the average transverse momentum (PxH^) P^^" 
portional to A^^^ [17] (and La oc A^^^ as well), one has that Axi oc A^^^ rather than A^^^. In 
what follows we consider two estimates for {p±w) [17]. The larger value, which comes from 
single nuclear rescattering of photoproduced dijets estimated in Ref. [18], is given by 



Considering that the initial states could not be explicitly identified, one assumes that (p^iy) 
is identical for quarks and gluons. In this case one has that when ~ 0.3 and A = 184, 
—dE/dz ~ 1.28 GeV/fm. We refer to this as "LQS" in the remainder of the discussion. The 
second estimate takes into account the difference between quarks and gluon interactions and has 
been derived in Ref. [17] considering the relation between the characteristic squared transverse 
momentum of the parton and the nucleon gluon distribution given by 



where xG{x) ~ 1—2 for the xi range of HERA-B. This lower estimate is referred to subsequently 
as "BDMPS". In this case one has that when a, ~ 0.3 and A = 184, -dEjdz ~ 0.12 GeV/fm 
and —dEg/dz ~ 0.28 GcV/fm. It is important to emphasize that a similar energy loss effect 
is expected for Drell-Yan production [4, 19]. In Ref. [19] the quark mean energy loss per unit 
length has been constrained to be —dEg/dz ~ 0.2 ± 0.15 GeV/fm considering a leading order 
analyzes of E866/NuSea and NA3 Drell-Yan data, which reasonably agrees with the BDMPS 
estimate. A comment is in order here. As pointed out in Ref. [4], the application of the 
BDMPS model a.t xp < may becomes problematic since Axi grows larger than xi, suggesting 
that the calculation may not be applicable for Axi > Xi. As it still is an open question, in a 
first approximation we will apply the model in the whole xp range. 

In order to analyze the quarkonium production we will consider the color evaporation model 
(GEM) [20]. In this model, the color charges of the cc produced arc randomized by the exchange 
of soft gluons, such that no information remains of the color configuration given by the preceding 
hard interactions. SU(3) algebra gives the probability 1/9 for the cc to be in a color singlet 
state and 8/9 to be in a color octet state. It is then assumed that all color singlet pairs with 
invariant mass below the threshold for open charm will form a charmonium state. The cross 
section for for the charmonium state i is then 



(pI^) ~ 0.658a,Ai/3GeV2 . 



(3) 



{pIw), = ^^PaxG{x,Q')La^ 0.07 as A'/' GeY' 
{pIw)9 = ^(pi^),^0.15a,AV3GeV^ 



(4) 




(5) 
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Figure 2: Effective exponent as a function of xp for cliarmonium production at HERA-B 
considering paladium and carbon targets. 

wfiere rricc is the invariant rnass of tlie cc pair, rric is tfie cliarm quark mass and 2mD = 3.74 GeV 
is the DD threshold. ^j'^J^ _ is the usual convolution of the perturbative QCD cross section 
with the parton density functions for the proton/nucleus, pi are nonperturbative universal 
factors which give the relative rates of producing the different charmonium states. Once pi has 
been determined for each state, e.g. ip, ip' or Xcj, the model successfully predicts the energy 
and momentum dependencies. We note that p^ includes both direct ip production and indirect 
production through radiative decays of the Xcj states and hadronic ip' decays. The pQCD cross 
section is taken in LO for simplicity, since we are just interested in ratios of cross sections. 
Also, the dependency on the overall factors pi/9 cancels out in the ratios. It is important to 
emphasize that although we will use the GEM to describe the quarkonium production, similar 
results are expected if the nonrelativistic QGD model (NRQGD) [21] is considered. 

In order to investigate the medium dependence of the quarkonium production cross section, 
we will follow the usual procedure used to describe the experimental data on nuclear effects 
in the hadronic quarkonium production [22], where the atomic mass number A dependence is 
parameterized by = CpN x Here and apjy are the particle production cross sections 
in proton-nucleus and proton-nucleon interactions, respectively. If the particle production is 
not modified by the presence of nuclear matter, then a = 1. A number of experiments have 
measured a less than linear A dependence for various processes of production, which indicates 
that the medium effects cannot be disregarded (See e.g. [5]). To estimate the modification 
of quarkonium production cross section due to the medium effects, we calculate the effective 
exponent a^xp), which is given by 

aM^L(^/!!^]/lr.A] . (6) 

[ \ axp I axp j I J 

Moreover, to obtain the available observable measured in the experiments, we also replace the 
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Figure 3: Effective exponent as a function of xp for cliarmonium production at HERA-B 
considering tungsten and carbon targets. 

nucleon by a liglit nucleus target (Carbon), calculating the ratio 

where Ai = C and A2 = Pd,W. Our results for the effective exponent aA2/Ai{xF) in J/^ 
production at HERA-B energy of ^/s = 41.6 GeV, with paladium and carbon targets are 
shown in Fig. 2 considering the EKS [Fig. 2 (a)] and DS [Fig. 2 (b)] parameterizations for 
the nuclear shadowing effects. In the dashed curves only the nuclear shadowing effects are 
taken into account. Energy loss effects are also included in the dot-dashed (BDMPS estimate) 
and full (LQS estimate) curves. The xp behavior of the effective exponent when considering 
the EKS and DS parameterization is similar. Whereas the shadowing alone produces a small 
enhancement for negative xp (antishadowing indeed) and a small suppression for higher positive 
values ofxp, the inclusion of energy loss leads to a large suppression for the heavy nuclear target, 
for negative values of Xi;'. This suppression can be understood in terms of the basic properties 
of the energy loss models. For a fixed value of negative xp, the corresponding value of x'l is 
smaller than the value of xi that enters in the parton distribution evaluation, since the projectil 
parton looses its momentum when traversing the nuclear target. This effect corresponds to a 
shift in the parton momentum to higher values of x in the nuclear case, which reduces the 
amount of partons in the initial state of partonic subprocesses as compared to the free nucleon 
scattering, since the parton distributions grow as x goes to zero due to the dynamical QCD 
evolution. 

In Fig. 3 we present our results for J/\E' production at same energy, in a tungsten and carbon 
targets, compared with preliminary HERA-B data for J/\I' production [7]. Concerning the EKS 
parameterization [Fig. 3 (a)], neither the enhancement due to the shadowing effects alone, nor 
the strong suppression for < due to energy loss are seen in the data, which indicate a 
small and a;i?-independent suppression. These results may indicate that the magnitude of the 
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Figure 4: Effective exponent as a function of for charmonium production at RHIC consid- 
ering gold and proton targets. 



antishadowing is overestimated in the EKS parameterization or that the energy loss is smaller 
that predicted by the LQS and BDMPS models. When the shadowing effects are taken into 
account via DS parameterization, the prediction get closer to the HERA-B data. This is due 
to the fact that the DS parameterization does not predict anti-shadowing behaviour of the sea 
quark and gluon nuclear distributions, as we can see in Fig.l. Again, when compared to the 
HERA-B data, the models BDMPS and LQS for the energy loss overestimate the suppression 
of the nuclear cross section for negative xp. 

Another possibility to constrain the magnitude of the medium effects is the study of the 
quarkonium production in proton(deuteron)-nucleus collisions at RHIC (See, e.g. Refs. [23, 
24]). In this case, due to the larger center of mass energy, the nuclear effects should be amplified. 
In Fig. 4 we present our estimates for the effective exponent for J/\I' production at \fs = 200 
GeV RHIC energy, with proton and gold targets. As we can see, the effective exponent is 
almost Xp independent in the positive xp range, with the EKS prediction being smaller due 
to the larger shadowing present in this parameterization. On the other hand, the energy loss 
leads to a very strong suppression for negative xp, larger than that predicted at HERA-B. It 
implies that the experimental analyzes of the quarkonium production is ideal to constrain the 
initial state energy loss effects in cold nuclear matter. 

One effect which we have disregarded in our analyzes was the nuclear absorption associated 
with the fact that the cc pair may interact with nucleons and be dissociated or absorbed before 
it can escape the target. This effect has been estimated in Ref. [25] considering different models 
for the quarkonium production and color singlet and color octet absorption. At HERA-B energy, 
the CEM and color singlet absorption in J/"^ production implies that a{xp = —0.5) ~ 0.97 
and a{xp > 0) ~ 1. On the other hand, if color octet absorption is considered one has 
a{xp = —0.5) ~ 0.96 and a{xp > 0) ~ 0.95. Consequently, if only the absorption effect 
is included in the calculations a reasonable description of the preliminary data is possible. 
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However, the inclusion of this effect in combination with shadowing and energy loss effects will 
implicate a larger suppression at negative xp, which is disfavoured by the data. Thus, the 
estimate of the nuclear effects in quarkonium production is still an open question. 

As a summary, in this paper we have studied the quarkonium production at HERA-B and 
RHIC. In particular, we have considered two distinct parameterizations for the shadowing and 
two estimates for the magnitude of initial energy loss effects and analyzed the Xp behavior of the 
effective expocnt. Our main emphasis was in the negative xp range which have been probed at 
HERA-B. We have verified that the inclusion of the energy loss strongly modify the behavior of 
a{xF) in this kinematical range. The comparison of our predictions with the preliminary data 
indicates that the combination of shadowing and energy loss effects, as described by the EKS or 
DS parameterizations and LQS or BDMPS models, is not able to describe the data. This may 
be related to the applicability of the BDMPS approach in cold matter and/or the overestimation 
of the antishadowing effect in the EKS parameterization. Another possibility is that the correct 
model for the initial state parton energy loss is one similar to the Gavin-Milana model [15], 
where Axi = eiXiA~' (e^ = 0.00412 and Cg = \eq), which implies a^xp) ~ 1 at negative xp. 
Another uncertainty present is the validity of the coUinear factorization in nuclear collisions. 
Some authors advocate that the coherence phenomena cannot be disregarded [26]. Since the 
suppression predicted in RHIC is larger, it would be desirable to have measurements in that 
region to constrain the different models and disentangle the different effects. To conclude, 
nuclear collisions remain a fascinating place to study different effects including the interplay of 
perturbative and non-perturbative QCD, and nuclear effects. 
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